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Abstract 

We show that at any location away from the spectral edge, the eigenvalues of the 
Gaussian unitary ensemble and its general f3 siblings converge to Sine/3, a translation 
invariant point process. This process has a geometric description in term of the Brow- 
nian carousel, a deterministic function of Brownian motion in the hyperbolic plane. 

The Brownian carousel, a description of the a continuum limit of random matrices, 
provides a convenient way to analyze the limiting point processes. We show that 
the gap probability of Sine/3 is continuous in the gap size and /3, and compute its 
asymptotics for large gaps. Moreover, the stochastic differential equation version of 
the Brownian carousel exhibits a phase transition at /3 = 2. 
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1 Introduction 



The Gaussian orthogonal and unitary ensembles are the most fundamental objects of 
study in random matrix theory. In the past decades, their eigenvalue distribution has 
shown to be important in several are as of probab ility combinatorics, number theory, op- 



erator algebras, even engineering (see lDeiftI (| 19991) for an overview). For dimension n, the 



ordered eigenvalues Ai < . . . < A„ G M have joint density 



e 



-/3EL=iAl/4 



where (3 = 1,2 for the Gaussian orthogonal and unitary ensembles, respectively. The 
above density makes sense for any /? > 0, and the point process is often called Coulomb 
gas in Gaussian potential at inverse temperature (3. The goal of this paper is to study its 
n ^ oo point process limit away from the spectral edge. 

The limit is described via a special case of the hyperbolic carousel. Let 

• 6 be a path in the hyperbolic plane 

• 2; be a point on the boundary of the hyperbolic plane, and 

• / : M+ be an integrable function. 

To these three objects, the hyperbolic carousel associates a multi-set of points on the real 
line defined via its counting function A^(A) taking values in Z U {—00, 00}. As time in- 
creases from to 00, the boundary point z is rotated about the center h{t) at angular speed 
A/(t). A^(A) is defined as the integer- valued total winding number of the point about the 
moving center of rotation. 

The Brownian carousel is defined as the hyperbolic carousel driven by hyperbolic 
Brownian motion h. See Section |2] for more details. 

In order to study the n ^ 00 limit of ^ we need to pick the center /i„ of the scaling 
window for each n. Then the scaling factor follows the Wigner semicircle law. Our main 
theorem gives necessary and sufficient conditions on /i„ to get a bulk-type limit. 

Theorem 1. For 13 > 0, let A„ denote the point process given by (H]), and let /i„ be a sequence so 
that n^/^{2y/n— +00. Then 

y/An- fil [An - fin) Sine/3, 

where Sine/3 is the discrete point process given by the Brownian carousel with parameters f{t) = 
(/5/4)e~^*/'^ and arbitrary z. 
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The convergence here is in law with respect to vague topology for the counting mea- 
sure of the point process. The limit and convergence for the special values = 1, 2, 4 



Mehta 



gwi) or 



under more re strictive scaling conditions has been well-studied, see 
Forrester! (|2008l) . The Brownian carousel description is novel even in these special cases. 
We note that the ensemble ^ may be generalized by replacing the J2k -^1 expo- 
nent by a similar sum involving a fixed function V of the eigenvalues. Assuming certain 
growth conditions on V the corresponding problem in the /3 = 2 c ase can be treated using 



ortho gonal polynomials and Riemann-Hilbert methods, see e.g. iDeiftI (|l999l) . 



Deift et al 



Together with the following theorem. Theorem [T] gives a complete characterization of 
the possible limiting processes for the ensembles ([T]). 



Theorem 2 (Ramirez, Rider, and Virag bOOTO ). For /3 > 0, let A„ denote the point process 
given by and let « sequence so that n/^{2y/n — fin) ^ a eR. Then 



n 



L/6 



(A„ - =^ Airjg + a 



Here Airy^ is defined as —1 times the poi nt pro cess of eigenvalues of the stochastic 
Airy operator, see Ramirez, Rider, and Virag (|2007l ) for more details. A straightforward 
diagonalization argument gives the following corollary, which is proved in Section |3l 

Corollary 3. As a ^ 00 we have 2y^(Airy^ +a) =^ Sine/?. 



The proo f of 



(l20o3) :)resent heuristics that the operators given by the tridiagonal matrices have a limit 



rotted il984) and 



heorem 1 is based on the tr i diagona l matr ix m odels introduced by 



Dumitriu and EdelmanI (|2002h 



SuttonI (|2005h and 



Edelman and Sutton 



whose eigenvalues give the Sine and Airy processes. Theorem |2] shows that this is in- 
deed the case at the spectral edge. The bulk case, however, is fundamentally different: 
there seems to be no natural limiting operator with the spectrum given by the Sine point 
process. Rather than taking a limit of the operator itself, we consider limits of discrete 
variants of the phase functions in the Sturm-Liouville theory. This connection is explored 
further in Section |53l where we describe how the Sine point process appears as a univer- 
sal limit for a large class of one-dimensional Schrodinger operators. 

The eigenvalue equation of a real tridiagonal matrix gives a three-term linear recur- 
sion for the eigenvectors. This becomes a two-term recursion for the ratios of consecu- 
tive entries, which then evolves by linear fractional transformations fixing the real line. 
So in our case, the evolution operators perform a time-inhomogeneous random process 
in PSL(2,M), the group of orientation-preserving isometries of the hyperbolic plane. To 
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Figure 1: The f3 = 1 stochastic sine equation as a function of A at three times 



get the Brownian carousel, we regularize th is evolution and take l i mits. An important 
tool is Proposition |23l (based on the results of IStroock and VaradhanI (|l979l )), which yields 



stochastic differential equation limits of Markov processes with heavy local oscillations. 

The Brownian carousel description gives a simple way to analyze the limiting point 
process. The hyperbolic angle of the rotating boundary point as measured from b{t) fol- 
lows the stochastic sine equation, a coupled one-parameter family of stochastic differen- 
tial equations 

dax = Xfdt + Re((e"*"^ - l)dZ), ax{0) = 0, (2) 

driven by a two-dimensional standard Brownian motion. For a single A, this reduces to 
the one-dimensional stochastic differential equation 

dax = \fdt + 2 sm{ax/2)dW, axiO) = 0, (3) 

which converges as t — oo to an integer multiple «a(oo) of 27r. A direct consequence of 
the definition of Sine/? is the following. 

Proposition 4. The number of points N{X) of the point process Sine/3 in [0, A] has the same 
distribution as ax{oo) / {211). 

Convergence to the solut ion of the coupled SDE s is the result formally a nnounced 
in the lecture by IVirag I (|2006l) . In independent work, iKillip and Stoiciul (|2006l) present a 
related b ut different desc ription of the limit processes in the setting of ci rcular ensern bles 
(see, e.g. iForresteil (|2008|) , Chapter 2 for discussion of these models and iKillipl (|2007|) for 
further related results). 

Proposition H] allows us to analyze the point process Sineg, for exa mple to deter mine 
the asymptotics of large gap probabil i ties. T his has been predi cted bv iDvsonI (119621) and 



WidomI |l996|) and for (3 = 1, 4 bvlTimbo et al.l (|198(3l) : there. 



proved for the cases /? = 2 by 
more refined asymptotics are presented; see also 

Theorem 5. For k > fixed and X ^ oo,we have 



Deift et al 



(|1997|) . 



P(# of points in [0, A] < A;) = exp ( - X'^{l3/6A + o(l))). 



This is shown in Section |Z3] for the case of a more general parameter /. Several similar 
asymptotic identities can be computed this way, and continuity properties can be studied. 
For the Sine/j processes we have 

Proposition 6. The probability distribution of N{X) is a continuous function ofX and (3. 

In contrast, the stochastic sine equation exhibits a phase transition at /5 = 2. 

Theorem 7. For any \ > Qwe have a.s. 

for all t large enough ax{t) > ax{oo) (4) 

if and only if (3 < 2. In particular, the probability of the event dU is not analytic at (3 = 2 as a 
function of f3. 

Deift (personal communication, 2007) asked whether this phase transition also ap- 
pears on the level of gap probabilities. This question remains open. 



2 The Brownian carousel and the stochastic sine equation 
2.1 Definitions 

In the Poincare disk model of the hyperbolic plane a boundary point can be described 
by an angle. The Brownian carousel ODE with parameters f{t) and zq describes the 
evolution of the lifted angle 7a (i) with e*'^^'^°^ = zq as it is rotated about the center B{t) at 
angular speed A/(t). Here B(t) is hyperbolic Brownian motion, that is the strong solution 
of the SDE 



B 



2^ 



dB = ^ '—^ dZ 

2 

driven by complex Brownian motion Z with standard real and imaginary parts. The 
speed of 7a, as measured in units of boundary harmonic measure from B, is \ f /(2tx). To 
change to an angle measured from 0, we need to divide by the Poisson kernel 



Poi(e^^, lu) = — Re 



27r e"i -w 27r |e*'^^ - w 
which yields the ODE 



2' 



Xf le*'^^ - 
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The most convenient way to define the winding number A^(A) of e''''^ about the moving 
center of rotation B{t) is to follow the corresponding angle. Let ax{t) denote the hyper- 
bolic angle determined by the points zq, B{t) and e*"^^*^*^. As we will check at the end of 
this section, Ito's formula shows that a satisfies the stochastic sine equation ([2]), i.e. 

dax = Xfdt + Re((e~^"^ - l)dZ), ax{0) = 0, (6) 

where dZ is simply complex white noise with standard real and imaginary parts. The 
name of the SDE comes from the fact that the last term equals 2sin(aA/2) Im(e~*"^/^(iZ). 
Since dW = Im(e~*"^/^(iZ) is 1-dimensional white noise, we get the SDE ^ for the single 
A marginals. 

Propostion|9]of the next section shows that 

^lima,(t) (7) 

exists for every A a.s. and for every Ai < A2 a.s. A^(Ai) < N{X2). Thus A^(A) can be defined 
as the unique random right-continuous function which agrees with ^ for every A a.s. 

To deduce ©, let T{w,z) denote the Mobius automorphism of the unit disk taking zq 
to 1 and taking w to 0. It is given by the formula 

^ S{w,z) z-w 
^{w^ z) = T^T r, S{w, z) = - — —. (8) 

S{W,Zo) 1 — wz 

Then a is defined as the continuous solution of 

a(0) = 0, e'°W =r(5t,e'^(*)). (9) 

The stochastic sine equation © follows from taking logarithms and applying Ito's for- 
mula. For the driving Brownian motion we get the explicit expression 

dZ = 2d2r{B, B)dB = ^— I dB. 

1 — \BY 1 — B 

Remark 8. By Ito's formula applied to the logarithm of ©, the noise term in ^ can be 
interpreted as the infinitesimal movement of the angle a under the difference of trans- 
formations dT = T{B + dB,'j)T(B,'j)^^. This infinitesimal Mobius transformation dT 
moves to T(B, B + dB) = d2T{B, B)dB, a standard complex Brownian motion incre- 
ment. Such a transformation dT changes the angle of any two points on the boundary by 
a Brownian increment with standard deviation proportional to their distance. This gives 
a more conceptual explanation of the noise term in ([6]). 
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2.2 Properties of the Brownian carousel 

Let Ll denote the set of absolutely integrable functions of IR+ which tend to at +00. 
Given a hyperbolic Brownian motion and a boundary point zq, the Brownian carousel 
associates a random counting function A^(A) to each f G L]. More generally, it is fruitful 
to study how A^(A) changes when the parameter / varies but the Brownian path remains 
fixed. In this case A can be absorbed in the parameter / so we will use the notation 
Nf = Nf{l), and «/ for the case A = 1. 

Proposition 9 (Properties of the Brownian carousel). We have 

(i) Of — Og has the same distribution as aj^g, 

(ii) a fit) is increasing in f, 

(Hi) [af{t)\ 2^ is nondecreasing in t when / > 0. Here [x\ 2^ = max (27rZ n {—00, x]). 

(iv) Nf = ^ linit^oo exists and is an integer a.s., 

(v) Nq = and Nf is increasing in f, 
ivi) E|iV^|< 2^||/||i, 

(vii) ENf = ^ f{x) dx, and 

(viii) Nf has exponential tails. For integers a,k > Owe have 



2na 



V{\Nf \ >ak)<2 

Proof. Claim [(1)] holds because «/ — solves 

dax = A(/ -g)dt + Re((e-^°^ - l)dZ*), ax{0) = 0, 
with dZ* = e~'^°'odZ. The standard coupling argument shows that the solution of the 



stochastic sine equation is monotone in the drift term, so we get (ii) 

Now assume that f > g = 0. Then by the above > 0!g{t) = 0. Claim (iii) follows 

by repeating this argument for the process after the hitting time oi2kn. 

Assume / > 0, and let F(t) = f{s) ds. Then — F is a continuous local martingale 

which is uniformly bounded below by — Thus it a.s. converges to a random limit. 
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So a also converges, but it can only converge to a location where the noise term vanishes; 
we get (iv)[ and (v) also follows from (ii) Now (vi) follows from 



Ea/(oo) - F(oo) < lim (Ea/(t) - Fit)) 



a;(0) - F(0) = 0, 



where the inequality is by Fatou's Lemma. By (iii) the function t ^ L'^/(^)J27r rionde- 
creasing, hence the above inequality implies that L'^/(^)J27r uniformly integrable, and 
so is Thus aj — F is a uniformly integrable martingale and so Eq;j(oo) = F(oo), as 



required for (vii) 



For general f e L], monotonicity |(ii)| gives 



< 



(10) 



where x"*" = max(x,0) and x = {—x)^. Now «-/- has the same distribution as —a/-. 
By the previous argument af+ and a/- are uniformly integrable. Hence aj — F is a uni- 



formly integrable martingale, and (iv) (vii) follow. Claim (v) also follows via (ii) We take 
positive and negative part of (|TOl) to get < af+, and < —a_f-, where the latter has 
the same distribution as «/- . Taking limits and expectations gives 



E|a/(oo)| = Ea;/(oo)^ + Ea/(oo) < Eq;/+(oo) + Eq;/-(oo) 



which gives (vi) 



[/a. Stopping 



Returning to / > 0, Markov's inequality implies that P(A^/ > a) < 
the process at time r when and if a hits 2Tcka we note that 

P {Nf >ik + l)a I Nf > ka, Tr) < 7^||M|i/a < ^||/||i/a, 
where is / shifted to the left by r. It follows that for integer k we have 

V{N^>ka)<{U\nVlo\' 
Ztx 

For general f E L\, we consider the positive and negative parts separately and use mono- 
tonicity [(v)] to get |(viii)[ □ 



Remark 10. The previous lemma shows that for a fixed f E L\ the random function 
A^(A) is a.s. finite, integer valued, monotone increasing with stationary increments. Thus 
A^(A) is the counting function of a translation invariant point process. Since f ^ aj and 
/ t-^ — tt-/ have the same distribution, the distribution of the point process is symmetric 
with respect to reflections. 
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Corollary 11. For any f E L\ the point process defined by N{X) is a.s. simple. 

Proof. The tail estimate of the lemma implies that the probability that there are two points 
or more in a fixed interval of length e is at most ce'^. Breaking the interval [0, 1] into pieces 
of length e, and using translation invariance, we see that the chance that there is a double 
point in [0, 1] is at most 2ce. Letting £ — > shows that a.s. there are no double points in 
[0, 1]. The claim now follows from translation invariance. □ 

Let M denote the space of probability distributions on Z with expectation. For /ii, /i2 G 
M let /i2) be the first Wasserstein distance, i.e. the infimum of E|Xi — X2I over all re- 
alizations where the joint distribution of {Xi, X2) has marginals fii and fi2. The topology 
induced by d is stronger than weak convergence of probability measures. Let jC{Nf) de- 
note the distribution of Nj. The following proposition is a stronger version of Proposition 
|6]in the introduction. 

Proposition 12. The map f t-^ J^{Nf) is Lipschitz-1 continuous in f: for f,g G Ll we have 
diCiNf),CiN,))<\\f-gh. 

Proof. Proposition |9] gives that Ng — Nj has the same distribution as Ng^f which implies 

d{CiNf), C{Ng)) <E\Ng- Nf\ = E\Ng^j\ <\\g- f\\,. □ 

2.3 Large gap probabilities 

Theorem 13. Let / : M+ ^ R+ satisfy f{t) < c/(l + t'^)for all t and \df \ < 00. Let k > 0. 
As X —>■ 00, for the point process given by the Brownian carousel with parameter f we have 

F{# of points in [0,A] < k) = exp {- XW\f\\l/8 + o{l))). (11) 

Lemma 14. Let Y be an adapted stochastic process with \Yt\ < m, and let X satisfy the SDE 
dX = YdB where Bt is a Brownian motion. Then for each a,t > Owe have 

P(X(t) - X(0) > a) < exp {-a^ /{2tm^)) . 

Proof. We may assume X(0) = 0. Then Xt = where r is the random time change 
r = Jq Y'^{s)ds. Since r < m^t the inequality now follows from 

P{Br >a) < exp (-aV(2r)) . □ 
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Proof of Theorem^3\ The event in (|TT|) is given in terms of the stochastic sine equation as 
limt_^oo ctxit) < 2k'K. We will give upper and lower bounds on its probability. 

Upper bound. By Proposition [9 | (iii) it is enough to give an upper bound on the probability 
that a stays less than x = 2{k + l)n. For < s < t we have 

P(a(t) <x\J^s) = P ^ 2 sm{a/2)dB > X fdt-x + a{s) J",^ . 

We may drop the a{s) from the right hand side and use Lemma [T4l with Y = —2 sin(Q;/2), 
m = 2, a = X{£ fdt-x/X) to get the upper bound 

( f* fdt - x/X)"^ 

P(a(t) < < exp(-AV(s,t)), r{s,t) = ^-^^ f^. 

6[t — S) 

Then, by just requiring a{t) < x for times e, 2e, . . . G [0, K] the probability that a stays less 
than X = 2{k + l)7r is bounded above by 

K/e K/e 

E JJP(a((A; + 1)£) <x\Tke) < exp | - ^ r(£A;, eA; + 

fc=0 fe=0 

A choice of e so that x/X = o{e) as A ^ oo yields the asymptotic Riemann sum 



y^r{ek,ek + e) = \ f f{t)dt + o{l). 



Letting K ^ oo provides the desired upper bound. 

Lower bound. Consider the solution a{t) of © with the same driving Brownian motion, 
but with initial condition a(0) = tt. Then a > a. For e < 7r/4, let As be the event that 

a{f) G (0, TT + e] for t G [0, s\. Then 

P(a(oo) < 27r) > P(A,) sup P(a(oo) < 27r | a(s) = y). 

S/e{0,7r+£) 

The sup is bounded below via Markov's inequality by 

Ti + e + Xr f{t)dt 
27r 

where the last inequality holds if s is set to be a large constant multiple of A. The event 
As is equivalent to i? = logtan(«/4) staying in the interval / = (— oo, logtan((7r + e)/A)\ 
where the evolution of R is given by Ito's formula as 

dR= - f cosh Rdt + -taxi\iRdt + dB, R{0) = 0. 
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Let /* = [—e, logtan((7r + e)/4:)], and consider a process R* so that (i) the noise terms of R 
and R* are the same and (ii) the drift term of R* at every time is greater than the spatial 
maximum over /* of the drift term of R. Let A* denote the event that for t E [0,s] we have 
Rl E I*. On this event R* > R, and therefore Ag also holds. With an appropriate choice 
of c we may set 



q{t) = (1/2 + ce)\f{t) + ce, 



dR* = q{t)dt + dB. 



(12) 



Let A* also denote the corresponding set of paths. Girsanov's theorem gives 

P(i?*eA:)=E l{BEA:)e^v[-\j\itfdt + j\{t)dB{t)^ . 
Integration by parts transforms the second integral: 

ns poo 

q{s)B{s)- B{t)dqit)>-C€\{f{s)+ \df\)-c£>-c'e{l + \) 
Jo Jo 

on the event B E A*^. Here we also used that / is bounded. The probability of this event, 
i.e. that Brownian motion stays in an interval of width ce, is at least exp(— c's/e^). In 



summary, (|T2|) is bounded below by 

exp {-ce{l + A) -cs/e^ - (1/8 + ce)X'^\\f\\l) . 
The choice s = cA, e = A~^/^ gives the desired lower bound. 



□ 



2.4 A phase transition at = 2 

The goal of this section is to prove Theorem [7] that at = 2 there is a phase transition in 
the behavior of the stochastic sine equation. 

As a converges to an integer multiple of 27r and it can never go below an integer 



multiple of 27i that it has passed (Proposition [9| |(iii)[ ), eventually it must converge either 
from above or from below. Theorem [7| says that a converges from above with probability 
1 if and only ii P < 2. Otherwise, it converges from below with positive probability. 

Proof. Case (3 < 2. It suffices to prove that if ax(to) E {2-nk — y, 2'Kk) with < y < tt then 
ax{t) leaves this interval a.s. in finite time. As a\{t + to) also evolves according to the 
stochastic sine equation with A' = Ae"'^*"/^, we may set to = and we are also free to set 
k = 1. Let B denote the event that the process ax started at ax{^) = x in (27r — y, 27r) will 
stay in this interval forever. It suffices to show that B has zero probability. 
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Consider i? = logtan(Q;A/4) andsety sothatlogtan((27r— y)/4) = 1. While « a ^ (0, 27r), 
Ito's formula gives the evolution of the process R: 



i?(0)=ro>l, dR = q{R,t)dt + dB, g(r, t) = ^ cosh r e"^*/^ + ^ tanhr 
Then B is the event that R{t) G (1, oo) for all t. On B we have 

q{R, t) > i tanh ^ > ^ - e"^^ > 1/4. 



(13) 



which gives R{t) > t/A - B{t) + tq from Set 

Q(t) = l ""'''f-' ds, (14) 
by the previous inequality, on the event B we have 

pt /"OO 

Q{t) > - / e-*/2+2B(s)^2ro^5 > _ / g-t/2+2B(s)-2ro^^ _ 

JO Jo 
where the integral M is a.s. finite. Let 

L{t) = R{t) - t/2 - B{t) - Q{t). (15) 

Then on B we have 

L'{t) = ^cosh{L{t)+t/2 + B{t) + Q{t))e-^'/^ (16) 

> ^exp[L{t) + Bit) + til/2 -P/A)-M]. (17) 

The equation follows from Ito's formula and the inequality uses coshr > e''/2. Multiply- 
ing (dZl by and integrating we get that on the event B 

e-m _ ^-m y(jf ^ ^^^/2 - /3/4)] ds. 

Jo 

with a random < C < oo. As the exponent is a Brownian motion with nonnegative 
drift, the limit of the integral on the right is a.s. infinite, thus the probability of B is 0. 

Case /5 < 2. It suffices to prove that for a large to if "aI^o) ^ (2vr — e, 2n) then axit) stays in 
the slightly larger interval (27r — 5, 2n) with positive probability. Choosing the values of 
e and 5 appropriately it suffices to show that if -R(O) > 2 and A is small enough then the 
event B that R G (1, oo) for t > 1 has positive probability. 

Recall the definition of Q and L from (|T4|) and (|T5|) . On the event i3 we have 

1/4 < tanh(/?) < 1/2, and - t/A < Qit) < 0. 
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Using this with (|T6|) and the fact that for r nonnegative cosh r < we get 



L'it) < ^exp[L(t)+5(t)+t(l/2-/?/4)]. (18) 

From ((T8|l we get 

e-^(o) _ e-^W < ^ Texp + s{l/2 - /5/4)] ds. 

Let M* denote the above integral for t = oo. Then M* is almost surely finite. Moreover, 
L{t) and thus R{t) remain finite if 

M*<2e-^(°V^- (19) 
From (dH) we get L'{t) > and > L(0) = 2 which gives 

Rit) > 2 + + t/2 + Q(t) > 2 + + t/4. 
So i?(t) stays above 1 if 

S(t)>-t/4-l foralH. (20) 

This has positive probability, so the conditional distribution of and M* given (|20)) is sup- 
ported on finite numbers. This means that the intersection of the events (|20l) and ((T9|) 
holds with positive probability for a sufficiently small choice of A, and it implies B. □ 



3 Breakdown of the proof of Theorem 1 



The goal of this section is to divide the proof of the main theorem into independent pieces, 
which in turn will be proved in the later sections. The proof presented here also serves as 
an outline of the later sections. 



Proof of Theorem^ Fix P > 0, and consider the n x n random tridiagonal matrix 

f Afo X(n-l)/3 ^ 
X{n-2)l3 
X(n-2)/3 J^2 



M(n) 



1 

7^ 



V 



(21) 



where the Xj/s ^rid J\fj entries are independent, J\fj has normal distribution with mean 
and variance 2, and Xjp has chi distribution with jP degrees of freedom. (For integer 
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values of its parameter, Xd is the length of a cZ-dimensional vector with independent stan- 
dard normal entries.) We let be the multi-set of eigenvalues of this matrix, which by 



Dumitriu and Edelman 



(|2002h has the desired distribution dT). 
First, we may assume that /i„ > 0; indeed, for a tridiagonal matrix, changing the sign 
of all diagonal elements changes the spectrum to its negative. In our case the diagonal 
elements have symmetric distributions, and by Remark [10] the limiting Sine/3 process is 
also symmetric. 



We set 



no = no(n) =n- ^1/4 - ^. 



The assumption n''^{2^fn — — > 00 implies 

no ^in ^ U, ^ i. 

So it suffices to show that 

if ^/i^/^ ^ as ?7. — i> 00, then 2?t,o''^(A„ — yU„) ^ Sine/3, (22) 

an equivalent version of the claim which makes computations nicer. Recall that the count- 
ing function A^(A) of a set of points in R is the number of points in (0, A] for A > or 
negative the number of points in (A, 0] for A < 0. 

Denote the counting function of the random multiset 2r;,|/^(A„ — yU^) by A^„(A), and that 
of Sine/3 by A^(A). Claim (|22l) follows if for every d>\ and (Ai, A2, . . . , A^) G M'^ we have 

(Ar„(Ai), Ar,(A2), . . . , iV„(Arf)) ^ (Ar(Ai), NiX^), N{Xa)) . 

The proof of this consists of several steps, these are verified in detail in the subsequent 
sections with the help of the Appendix. 

Consider the one-parameter family of SDEs defining the Sine/3 process: 

dax = A^e-^*/^dt + Re((e-*"^ -l)dZ), (23) 

where Z is complex Brownian motion on [0, 00) with standard real and imaginary parts. 
The time-change t — | log(l — t) transforms (|23l) to 

2^/(3{l-t) dax = XP'^'dt + 2V2Re((e-^°^ - l)dW), (24) 

where Wt is complex Brownian motion for t E [0, 1) with standard real and imaginary 
parts. Proposition |9] of Section |2] shows that the counting function A^(A) of the process 
Sine/3 can be represented as the right-continuous version of (27r)"Mim(^oo c^aI^)/ a limit 
which exists for every A G M a.s. This gives 
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Step 1. For every A G M, a.s. we have 27rA^(A) = lime^o+ c^xi^ ~ 

The eigenvalue equation for a tridiagonal matrix gives a three-term recursion for the 
eigenvector entries. This can be solved for any value of A, but the boundary condition 
given by the last equation is only satisfied for eigenvalues. 

The ratios of consecutive eigenvector entries a evolve via transformations of the 
form r i-^ a — b/r,b > 0. These transformations are isometries of the Poincaree half plane 
model of the hyperbolic plane. The hyperbolic framework is introduced in Section |4lT] for 
the study of these recursions. In particular, a moves on the boundary of the hyperbolic 
plane which can be represented as a circle, eigenvalues can be counted by tracking the 
winding number of a as a function of A. The rough phase function a (introduced in 
Section |42|) transforms r^ A to an angle through 2 arctan(r^ a)- Taking always the appro- 
priate inverse of tan we get a continuous function of A taking values in M, the universal 
cover of the circle. 

Our goal is to take limits of the evolution of cpe^x. Since it has fast oscillations first it 
needs to be regularized. In order to remove the oscillations, we follow a shifted version of 
the hyperbolic angle of a around the fixed point of a simplified version of the transfor- 
mation r a — b/r. As we will see later, the important part of the evolution takes place 
in the interval < f < [rioj , which is exactly when this transformation is a hyperbolic 
rotation. 

The precise regularization is done in Sections |4.3t there we introduce the (regularized) 
phase function cpe^x and target phase function ipf^ with parameters < i < [uq] and 
A G M. These correspond to solving the eigenvalue equations starting from the two ends, 
1 and n. As n ^ oo, these two parts will require completely different treatment, so it is 
natural to break the evolution into two parts this way. Proposition [18] shows how we can 
count the eigenvalues using the zeroes of these phase functions mod 27r; this is a discrete 
analogue of the Sturm-Liouville oscillation theory. 

Let denote the number of elements of A. 

Step 2. For i = 1, 2, ... , [^oj, the function a monotone increasing, and is independent of 
(pfy For any A < A' and I < i < [noj almost surely we have 

iVn(A') - iV„(A) = # {iipe,x - iflx. - ^lx^ n 27rZ) . (25) 

Since (p counts all eigenvalues below a certain level, its regularized version ip will 
encode all the fluctuations in the number of such eigenvalues. So the continuum limit of 
V9 is expected to have large oscillations as its time-parameter converges to oo. In order to 
deal with this problem, we introduce the relative phase function a = — V^£,o- This is 
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related to the number of eigenvalues in an interval, so it is expected that its scaling limit 
will have nice behavior at +00. Note that ai^x has the same sign as A by Step|2l Let 

mi = [no{l - e)\, 1712 = [n - - k (/i^^' V 1)J, 

where the constants £,k > will be specified later in a way that the chain of inequalities 
< mi < m2 holds. 

Next we will describe the limiting behavior of ^pi^x and a£^x when i is in the intervals 
[0, mi] and [mi, 1712], respectively. This is the content of the next three steps. Section |5J] 
studies the behavior of the relative phase function on [0,mi]. In Corollary |27| we will 
prove that ae^x converges to the SDE (|24|) in this stretch. 

Step 3. For every < £ < 1 

amux==^ ax{l - e), as n 00 (26) 

in the sense of finite dimensional distributions for A. 

Proposition |28] of Section 16.11 shows that ai^x does not change much in the second 
stretch, if it is already close to mod 27r at the beginning of the stretch. 

Step 4. There exists constants cq, ci depending only on A and f3 such that if k = Kn > cq, X < |A| 
then 

E [\{am„x - «™2,a)| a 1] < Ci(Edist(a^,,A, 27rZ) + e'^' + V 1) + k''), (27) 

In Proposition |33] of Section [63l we show that (fmzfl becomes uniform mod 2-k. 
Step 5. If K 00 and riQ^KdMl"^ V 1) — * then 

{Vm2,o}27T =^ Uniform[0,27r], 

where {x}^^ = mmk(^z,k<xix - 27rfc). 

Finally, in Lemma |34] of Section [6T4l we show that nothing interesting happens after m2. 
Step 6. For every fixed k > and A G M 

k®2,A-<^m2,o| -^0. asn^cx). 

In a metric space, if lim„^oo Xn,k = Xk for every k and also lim/j^oo Xk = x then we can 
find a subsequence n{k) 00 for which lim„^oo a;„,„(fc) = x. This simple fact, together 
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with the previous steps, allows us to choose sequences £ = — > 0, k = k„ ^ oo in a way 
that the following limits hold simultaneously: 

iam„x,,t = l,...,d) ^ 2n{Ni\i),i = l,...,d), (28) 
Wm,,o}2^ ^ Uniform[0,27r] (29) 
Ka-<,o| ^ 0, ^ = l,...,d (30) 

Since dist(-, 27rZ) is a bounded continuous function, (|28l) implies that the right hand side 
of (|27|) vanishes in the limit and so 

ttmi.A, - ttma.A, 0, i = 1, . . . , (31) 

By (|28ll and dST) the completion of the proof only requires the following last step. Let {x)^^ 
denote the element of 27rZ in [x — vr, x + tt). 

Step 7. For i = 1,. . . ,d and X = XiWe have lim„_oo P (27rA^„(A) = (ttm2,A)27r) = ^■ 

We conclude by the proof of Step |Zl We suggest skipping it at the first reading, as it is 
the most technical part of this outline. We include it here because it uses too much of the 
notation and assumptions of the preceding discussion. 

We will assume A > 0, the other case follows similarly. Then < {am2,x)2n ^ 27rZ, for 
any x G M we have 

= 27r#((x,x + (a^,,A)2j n 27rZ) 
Using this with x = y^^^.A - ^%,,x - («m2,A)2^ we get 



)2^ = 271 # (<^m2,0 - 2,A + "m2,A - (am2,A)2^ , V'm2,A " 2,a] ^ 27rZ) . (32) 



The symmetric difference between the intervals in (|32l) and (|25l) is an interval J with end- 
points v?m2,o - ^^2,0 and v?m2,o - V5^2,A + "m2,A " («m2,A)2^- So it suffices to show 

lim P((Jn 27rZ) = 0) = 1. (33) 

n— >oo 

We will show that the length of J converges to while one of its endpoints becomes 
uniformly distributed mod 2%. First, 



1^1 < |am2,A - (ttm2,A)2^| + |</'m2,A " ^m2,o| ^ 



where the convergence of the first term follows from (|28l IST]) as a.m2,x converges to an 
element of 27rZ; the convergence of the second term is (|30l) . Also, since (pm2,o arid v^®^ o are 
independent, from (|29|) we have {fm2,o — fm2fi}2n Uniform[0, 27r]. Equation (|33l) . Step 
[7| and the theorem follows. □ 
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Proof of Corollary\3\ Note that weak convergence of point processes is metrizable. Let — > 
oo. For every i, we can find Ui > i so that the point process 

A* = 2y^nf (A„, - + aim-'/') 

is 1/i-close to 2y/a~(Airy^ +aj) by Theorem|2] By Theorem[T} A* converges to Sine/3. □ 

4 The hyperbolic description of the phase evolution 
4.1 The hyperbolic point of view 

The eigenvector equation for a tridiagonal matrix gives a three-term recursion in which 
each step is of the form m^+i = bug — au^^i, in our case with a > 0. Let PSL(2, R) denote 
the group of linear fractional transformations preserving the upper half plane HI and its 
orientation. Then = ug^i/ui evolves by elements of PSL(2, M) of the form r ^ h — a/r. 

We will think of EI as the Poincare half -plane model for the hyperbolic plane; it is 
equivalent to the Poincare disk model U via the bijection 

i + z 

which is also a bijection of the boundaries. Thus PSL(2, M) acts naturally on U, the closed 
unit disk. As r moves on the boundary dW = M U {00}, its image under U will move along 

dV. 

In order to follow the number of times this image circles U, we would like to extend 
the action of PSL(2,M) from dU to its universal cover, W = M, where we use prime to 
distinguish this from dM. This action is uniquely determined up to shifts by 27r, but here 
we have a choice. For each choice, we get an element of a larger group UPSL(2, M) defined 
via its action on M'. UPSL(2,M) still acts on M and U just like PSL(2,M), and for T G 
UPSL(2, M) the three actions are denoted by 

M^M: z.T, U ^ U : ^ ^ ZoT, R' ^R' : zh^ z*T. 

We note in passing that the topological group UPSL(2,M) is the universal cover of the 
hyperbolic motion group PSL(2, M), and PSL(2, R) is a quotient of UPSL(2, R) by the in- 
finite cyclic normal subgroup generated by the 27r-shift on R'. For every T G UPSL(2, R) 
the function x t— > a;*T is strictly increasing, analytic and quasiperiodic, i.e. {x + 27r)*T = 
x*T + 27r. 
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Given an element T G UPSL(2, M), x,y E M', we define the angular shift 

asliffi/(T, X, y) = {y*T - x*T) - {y - x) 

i.e. the amount the signed distance of x, y changed over the transformation T. This only 
depends on the image of T in PSL(2, R) and the images v = e", w = e*^ G dU of x, y under 
the covering map. This allows us to define ash(T, v, w); more concretely, 

ash{T,v,w) = asli]R/(T,a;,y) = Arg^^^^^^{waT/voT) - Arg^^^^^^{w/v), 

where the last equality has self-evident notation and is straightforward to check. Note 
also that the above formula defines a.sh.{T,v,w) for T G PSL(2,]R), v,w G dV as well. 
For explicit computations, we will rely on the following fact, whose proof is given in 
Appendix lA.il 

Fact 15 (Angular Shift Identity). Let T G PSL(2,R) be a Mobius transformation and v,w e 
dVjleta = O.T-\ Then 



ash(T, V, w) = 2Arg 



,W — (7)V 



w[v — a 



, 1 — aw 
2Arg( - 



(34) 



Next, we specify generators for UPSL(2,M). Let Q(a) denote the rotation by a in U 
about 0, more precisely, the shift by a on M': 



(p*Q{a) = + a 



(35) 



For a, 6 G M let A(a, h) be the affine map z t-> a(2; + 6) in H. If a > then this is in PSL(2, M), 
it fixes the oo in 9EI and —1 in dV. We specify the action of A on M' by making it fix tt G M'. 
Then we have 

^*A(a, 6) = + ash(A(a, 6), -1, e^'^). (36) 
The following lemma estimates the angular shift. The proof is given in Appendix lA.il 
Lemma 16. Suppose that for a T G UPSL(2, R) we have [i + z).T = i with \z\ < 1/3. Then 



ash(T, V, w) 



Re [{w - v) i^-z - ^^^±f±^ 

— Re [{w — v)z] + 62 
£1, 



+ 63 



(37) 



where for d = 1, 2, 3 and an absolute constant c we have 

kd| < c\w - v\\z\'^ < 2c\z\'^, 
Ifv = —1 then the previous bounds hold even in the case \z\ > 1/3. 



(38) 
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4.2 Phase evolution equations 

The eigenvalue equation of a tridiagonal matrix can be solved recursively. The goal of 
this section is to analyze this recursion in terms of phase functions. 

We conjugate the matrix M = M{n) in (|2T|) by a diagonal matrix D with 



Du = Din) a 



n 



X(n-£)/3 



where 



^3 



\/n- i ~ 1/2. 



.1 

We get the tridiagonal matrix = D~^MD given by 

si Xi si + Yi 



1 

7^ 



M2 



\ 



■ J 



S2 



\ 



(39) 



Then M and M have the same eigenvalues, but M has the pro perty that the eigen 
value equations are independent. (A similar conjugation appears in 
(|2007|) .) The moments of the independent random variables 



Edelman and Sutton 



X, 



X(n-j-l)/3 



— S 



i+1 



J' 



< J < n - 1 



are explicitly computable via the moment generating functions for the F distribution. 

Our proof is valid for any choice of independent real-valued random variables Xj, Yj 
satisfying the following asymptotic moment conditions. Xj and Yj may also depend on 
n, in which case the implicit error terms are assumed to be uniform in n. 



(40) 



Let M£ = A (1 < ^ < n) be a non- trivial solution of the first n — 1 components of the 
eigenvalue equation with a given spectral parameter A, i.e. 



moment 




2nd 


4th 






2/P + Oiin-j)-') 


0{1) 



SiUi + Xi Ui+i + (Yg + Si)ui+2 = 

Uo = 0, Ui = 1. 



< £ < n - 2 



(41) 



Then with = a = ug^i/ug we have 



re+i 



A Xf 



+ 1 + - 



Y 



< £< n-2 



(42) 
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This also holds for = or = oo; in fact, the initial value of the recursion is tq = oo. If 



we set Yn-i = and define via the £ = n — 1 case of (|42|), then A is an eigenvalue if and 
only if r„ = 0. 

We will use the point of view and notation introduced in Section |4H Namely, r takes 
values in dM = M U {oo}, the boundary of the hyperbolic plane. Moreover, the evolution 
of r can be lifted to the universal cover of dM. The extra information there allows us to 
count eigenvalues, as the following proposition shows. The proposition also summarizes 
the evolution of r and its lifting G M'. We note that this is just a discrete analogue of the 
Sturm-Li ouville oscillation theory suitable for our purposes; such analogues are available 
in the literature. Although we state this proposition in our setting, a trivial modification 
holds for the eigenvalues of general tridiagonal matrices with positive off-diagonal terms. 

Proposition 17 (Wild phase function). There exist functions (p^(p® : {0, 1, . . . , n} x R ^ M 

satisfying the following: 

(i) re,A.'U = e'*^'^, 

(ii) <^0,A = TT, (^®^ = 0. 

(Hi) For each < i < n, (pi^t^ is analytic and strictly increasing in A. For < £ < n, (p®^ is 
analytic and strictly decreasing in A. 

iiv) For any < £ < n, A is an eigenvalue of M if and only ifffe^A — <^®a ^ ^vrZ. 

Proof. We consider the following elements of the universal cover UPSL(2,M) of the hy- 
perbolic motion group PSL(2, M): 

Q(7r), W, = Aiil + Y,/s,)-\-Xj/sj) 0<j<n-l (43) 

where Q corresponds to a rotation in the model U, and A corresponds to an affine map in 
the model H, as defined in (|35ti36[) . With this notation, the evolution (|42]) of r becomes 

R,,A = Q(7r)A(l,A/s,)W,, ^^^^ 

for < i < n — 1, and A is an eigenvalue if and only if oo.Rq^a ■ ■ • Rn-i,A = 0. Multiplying 
this by (R£,A ■ ■ ■ Rn-i,A)"^ for some < i < n and then moving to the universal cover M' 
of dM gives the equivalent characterization <^£_a = mod 2%, where 

'p£,A = 7I"*Ro,A ■ ■ ■ R^-1,A, 'ff^A = 0*-f^n~l,A ' ' ' -^^^A' i'^^) 

which is exactly (iv). Claims (i)-(ii) follow from the definition. 
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As v9o,A = TT, one readily checks that v^i a is strictly increasing. Since {ip, A) i-^ (/'*R^,a 
are nondecreastng analytic functions in both parameters and so are their compositions, 
the statement of claim (iii) for a now follows. The same proof works for (^®. □ 

Motivated by part (iv) of the proposition, we call ip® the target phase function. 



4.3 Slowly varying phase evolution for a scaling window 

For scaling, we set 

s(r) = s(")(r) = ^l-T- 1/2" 

so that we have = s{i/n)y/n. Making s depend on n via the ^/2n term helps make the 
upcoming formulas exact rather than only asymptotic. 

The phase function (p£ introduced in the previous section exhibits fast oscillations in i. 
In this section we will extract a slowly moving component of the phase evolution whose 
limiting behavior can be identified. The oscillations of are caused by the macroscopic 
term Q(7r) A(l, A/s^) of the evolution operator R^,a- The recursion (|44|) has different be- 
havior depending on whether this macroscopic part is a rotation or not. As we will see 
later, the continuum limit process comes from the stretch < i < uq where it is a rotation 
(this is because the corresponding eigenvectors will be localized there). The eigenvalues A 
of interest will be near the scaling window fin, so we define the main part of the evolution 
operator as the macroscopic part of R^,;^^, that is 

J, = Q(7r) A(l,/iJ.,) = Q(7r)A(l, J"; ). (46) 

This is a rotation if it has a fixed point pi in the open upper half plane H, the fixed point 
equation p^.^i = pe turns into 



Note that v^ira is the relative location of the scaling window in the Wigner semicircle 
supported on [—1,1]. Since s(r) is decreasing, we have that pi e M for r < no/n, where 
s{nQ/n) = pn/ V4n. This explains the choice of the parameter tiq. 

Thus pi = p{no/n, £/n), where p(ri, T2) is the solution in the closed upper half plane of 



/-24^P+1 = 0, i.e. P(-i,-2) = 4^ + Jl-4^. (48) 

s{T2) S{T2) V s{T2y 



More specifically. 



/ Pn/'^ ■ n-o — i 
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Because of our choice of scaling window and the density in the Wigner semicircle law 
it is natural to choose the scaling (|22|) by setting 

^ = ^ln + TT^- (50) 



We recycle the notation ui^\, ri^x, (pi^x-, <^®a f*^^ the quantities m^^a, r^^A, '^i,A-, ^f^A- separate 
Ji from the evolution operator R to get: 

R„ = J,L,,W„ L,, = a(i.5-^^^^). (51) 

Note that a and become infinitesimal in the n — > oo limit while Ji does not. Ji is a 
hyperbolic rotation, differentiating z t— > z.Ji at z = pi shows the angle to be — 2Arg(p^) G 
[-7r,0]. Let 

T, = A(Im(p,)-\-Re(p,)) 
correspond to the affine map sending pi E M to i E M, then we may write 

J, = Q(-2Arg(p,))T", 

where = B^^AB. Rather than following (p itself, it will be more convenient to follow 
a version which is shifted so that the fixed point pi of the rough evolution is shifted to i. 
Moreover, in order to follow a slowly changing angle, we remove the cumulative effect 
of the macroscopic rotations J^. Essentially, we study the "difference" between the phase 
evolution of the random recursion and the version with the noise and A terms removed. 
The quantity to follow is 

V^^A = 'p£,x*TiQ£-i, (52) 

where 

Q, = Q(2Arg(po)) . . . Q(2Arg(p,)), -1 < £ < no- 
Acting on U, is simply a rotation about 0, more precisely a multiplication by 

r]e = pIpI . . . pI (53) 

From (|45|) and (|52l) we get that cp evolves by the one-step operator 

(T,Q,_i)-1R,,a(T,+iQ,) = (T,-iL,W,T,+i)Q^ := (S,,a)^^ 

We keep this "conjugated" notation because S^^a corresponds to an affine transformation. 
For i < uq we define the corresponding target phase function 

= <^®,.T,Q,_i. (54) 

The following summarizes our findings and translates the results of Proposition [TTIto this 
setting. Here and in the sequel we use the difference notation Ax£ = x^+i — X£. 
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Proposition 18 (Slowly varying phase function). The functions v?, V'® : {0, 1, • • • , L'^oJ} x 
M ^ R satisfy the following for every < £ < uq: 

(ii) (fii^x ^nd —(pfx i^^s analytic and strictly increasing in X, and are also independent. 
(Hi) With S£,A = T^^L^W^T^+i, we have Aipe,x = ash(S£,A, -1, e''f''^f]i). 

(iv) ^e,x = ^iMQili^i^- 

(v) For any X < X' we have a.s. Nn^y — Nn,x = # ((v'^.a — V^^A' V^^,a' — vfx'] ^ SttZ). 
The form 

breaks S into a deterministic A-dependent part and a random part that does not depend 
on A. Let a = V^£,a*(L^,a)^'^^^ be the intermediate phase between these two steps. Note 
that (p}q = ^e,o> arid 

(L,a)^^ = Afl, ^ ^ =)=Afl, , ^ ) (55) 

The relative phase functions 

ae,x = ^e,x - ^e,o, a^x = vlx - Vlo 
are the main tools for counting eigenvalues in intervals. 

Proposition 19 (Relative phase function). The function a : {0, 1, ... , [rioj } x M ^ M satisfies 

(i) C(o,x = 0, a£fl = and for each i > 0, a^^x is analytic and strictly increasing in X. 

(ii) Aae^x = asli((LA,„-f)^% -1, e^^^^^r/^) + ash(Sf,o, e^^^*^?]^, e^'^^-oy^^) 

= ash((L,,A)^^ -1, e^^^-^r/,)+ash(S,,o, e'^^l^f],, e^^^.^r/,)+ash(S,,o, e'^^-^r/,, e^^*-"//,) 

(Hi) For each i and X > Owe have [«£,Aj27r ^ I_'^^+i,aJ ^ L'^^+i,Aj27r- 
Proof (i)-(ii) are direct consequences of Proposition [181 To check (iii), we note 

cxe,x = fe,x - fe,o 
<^lx = Ve,x*{'^e,x)'^'^' - y^e,o 
«m,A = ¥^^,a*(L,,a)''^'^^(S,,o)'^^-¥^£,o*(S^,o)'^^ 

Since the map a and its conjugates are monotone in A, we get a < a- Since (S^^o)*^^ 
is the lifting of a Mobius transformation, it is monotone and 27r-quasiperiodic, whence 

NJa^ = bm.AJa^- □ 
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Remark 20 (Translation to the original matrix). Let denote the solution of the discrete 
eigenvalue equation for the original matrix (|2T]) . It is given in terms of the diagonal matrix 
D defined in the beginning of the section and the solution u of recursion (|4T|) as wg = 
{Du)i. The ratios of the consecutive entries of this vector are 

wi+i {Du)i+i Ui+i Di+i,i+i X{n-e-i)i3 
We {L>u)e Ui Di^i VPSe-i 

If i < no then we may further rewrite this using zi as 

pe = ^(^^ziM ((^,.U-^)r/,_i Im(p,) + Re(p,)) . 

VP se+i 

4.4 The discrete carousel 

Corollary |27| in Section \52\ shows that the appropriate limit of the relative phase function 
a£^x is the stochastic sine equation. In this section we bring the discrete evolution equa- 
tions in the form that it becomes clear that their limit should be the Brownian carousel. 

By (|44l) the evolution of is governed by a certain discrete process G^^a in the hyper- 
bolic automorphism group UPSL(2, R): 

re = ro.Gtx = '^o-R-ca ■ ■ ■ R.£-i,a- 

This process has rough jumps, but it is a smooth function of the parameter A. It is therefore 
natural to expect that the evolution of the automorphism a^7o ^^^^ have a continuous 
scaling limit. In the following, we will rewrite this expression in a form indicating the 
desired scaling limit. 

By (|52|) the evolution operator G^ a of satisfies G^ a = G^^aT^Q^^i, and therefore 
G^ A^7o ^ x^eo- evolution of (pe,\ is given by 

= V5o,A*G£,A = T<'*Ge,x, (56) 

where with = B~^AB we have 

G^^A = Yq^aXq Yi^aXi ■ ■ ■ Y^_i^aX^-i 

- J^o.A I i,A ^2,\ ■ ■ ■ ^e- 
Gi = Gift = XoXi---X£_i, (57) 

and we used the temporary notation a = i{1^e,x)^^)^', = (Scfi)^'. By definition, 

ai,x = ^eMQ{-^e,o) = tt^G^^a Q{-^e,o)- (58) 
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We introduce the notation 

7£,A := vt.CaG,-! = 7r*Yo,A Y^f vff ■ ■ ■ vfj^), (59) 
Be := OoG7^ G U. (60) 

With T denoting the Mobius transformation defined in ((8]), we claim that 

T{Bi, z) = ZoGt Q(-V5£,o) 

with the choice of zq = —1. This follows from the fact that T{Be, Bi) = by definition and 
T{Bi, -1) = 1 by (|5El5Zl»- Hence ^ becomes 



which is the same form as equation (HJ) relating the stochastic sine equation to the Brown- 
ian carousel ODE. 

Remark 21 (Heuristics). Note that is approximately an infinitesimal noise element in 
PSL(2, R). acting on U moves infinitesimally in a random direction. This direction 
is not necessarily isotropic, but the conjugation by the macroscopic rotation makes the 
composition of consecutive X^'s move to an approximately isotropic random direction. 
Thus Bi in (|60|) approximates hyperbolic Brownian motion started at run at a time- 
dependent speed. Similarly, the a are infinitesimal parallel translations, but because of 
the conjugation by the macroscopic rotations Q^, their composition approximates rotation 
about 0. Thus 7^^a in (|59l) approximately evolves by rotations about Be. This is exactly 
how the Brownian carousel evolves, giving a conceptual explanation of our results. This 
suggests an alternative way to prove our results via the Brownian carousel ODE (O. 



5 The stochastic sine equation as a limit 

This section describes the stochastic differential equation limit of the phase function on 
the first stretch [0, no(l— £:)]. In the limit, this stretch completely determines the eigenvalue 
behavior; this will be proved in Section [6l 

5.1 Single-step asymptotics 

Let JF£ denote the cr-field generated by the random variables Xo,Xi, . . . , X^.i, andFo? ^i? • • • ,Ye-i- 
Let E£[- ] denote conditional expectation with respect to JF^. By definition, the random 
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variables <^£^a, '^i-, di are measurable with respect to JF^. Moreover, for fixed A, both a 
and i-pi^x are Markov chains adapted to J^i. 

Throughout this and the subsequent sections we assume that |A| is bounded by a con- 
stant A. By default the notation 0{x) will refer to a deterministic quantity whose absolute 
value is bounded by c\x\, where c depends only on /3 and A. As £ varies k will denote 
no - i. 

This section presents the asymptotics for the moments of step A(/9£ a := — ^f.,x- 

Recall from Section |43l that ^ moves on the interval [0,r2o]. The continuum limit of i^^^x 
will live on the time interval [0, 1] so we introduce 

t=-Gp,l. 

no 

We also introduce the rescaling of s(t)^ on this stretch: 

..^.2 _ s{tno/nf - s{no/nf ^^^^ 

with s > 0. This actually simplifies to 



s{t) = VT^ = ^/kj^o (62) 

in our case. As we will see later, the scaling limit of the evolution of the relative phase 
function will depend on s and the scaling parameters through s. The fact that this func- 
tion only depends on t explains why the point process limits do not depend on the choice 
of the scaling window in Theorem [U In Section 15.31 we provide a more detailed discus- 
sion and further implications. We will keep the notation s (instead of writing a/1 — t) to 
facilitate the treatment of a more general model discussed there. 

Proposition |18|(iii71 expresses the difference Acpi^x := ipe+i,x — fi,x via the angular shift 



of A arid (L^ a)^*- Lemma [T6l in turn, writes the angular shift in terms of the pre-image 
ofi = a/— 1. In the present case 

Ze,x = ^.Sll - t = 2.T7_^\(L,,aW,)-^T, - z = ve,x + V,, (63) 

where 

^ , P£+i - Pi Xi + pi+iYi 

Vi,X = TT-T + — f , Ve = — . (64) 

2nos{t) imp^ \/no s[t) 

The random variable Vi is measurable with respect to but independent of J^^. 

By Taylor expansion we have the following estimates for the deterministic part of Z^ ^: 
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where we abbreviate p{t) = p{no/n,tno/n) = pe, see (|48|) . The behavior of the random 
term is governed by 



EV? 



0(n - i)-'/'k~"^' 

q{t) + 0{n ~ (.)-y^k~'''\ E\Vi\'^ = 0{k"^'^ 



E\Vl\ = ^^p{t) + 0{n-e)-'k 



J_ 

no 



d = 3A 



where 



Pit) 



4nn 



2(1 + p2 



(66) 



(67) 



ps^ pk' ' ' 

Here the error terms come from the moment asymptotics (|40l), the size of s, and from the 
bounds 



p,+,-p, = Oin~i)-y'k-'/' 



d pe+i - pe 

T+P 

at riQ 



0{n - ty^i^k'^^ 



Proposition 22 (Single-step asymptotics for v?£,a)- for i < uq with t = i/riQ and k = uq — 

we have 



E [A(pi^xAipi^x'\ve,x = x,(pi^x' = y] 



1 



bx{t) + —osci + 0{k''^) = 0{k 



no no 

—a{t, X, y) + —0SC2 + 0{k 

no no 



3/2 X 



(68) 
(69) 



E^ |A^. 



0{k 



-d/2-s 



d = 2,3, 



where 



2s Imp ' 2PP ' 
The oscillatory terms are 

osci = Re {{-vx - iq/2)e~'''r]i) + Re (ie"^"'?7f g) /4, 

0SC2 = pRe(e-"r/£ + e-*^r7£)/2 + Re(g(e-*> + e-^^r/£ + e-^(^+^)r/f))/2. 



(70) 



(71) 



Proof. By Proposition ll8|(iiT) the difference A(/9£ a can be written 

Aipi^x = ash(S£,A,-l,^^?) 



Re 



-{l + zr])Z- 



4 



+ C(Z3 



(72) 



ReZ -' ' ^ ' ^'^''^•^^ 



4 



+ 77 terms + 



where we used Z = Zi^x> V = Ve ^rid z = exp{iipc^x)- The estimate (|72)) is from the quadratic 
expansion (|37|) of the angular shift in Lemma [161 Note that since the second argument of 
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ash is —1, we do not need an upper bound on \Z\. We take expectations, the error term 
becomes 

0{E\Z\') = 0{\ve,x\' + E\Ve\') = Oik'''). 

By dllUlll) we may replace the EZ , E|Zp and EZ^ terms by vx{t) , p{t) and q(t) while 
picking up an error term of 0{k~'^). Significant contributions come only from the non- 
random terms v^^x of Z and the expectation of V^. We are then left with oscillatory terms 
with T], error terms, and the main term 

-Ref^,A + ImEV//4 = — {-Revx + lr^vq) + 0{k''^) 



no 



1 

no 



A Rep' ^ \m{p' 



2s Imp 2/?s2 



+ 0{k 



The error terms come from the moment bounds (|40|) of X, Y and from the discrete ap- 
proximation of the derivative Re p'; their exact order is readily computed. This gives (|68|) 
with (|7T]) . The C(A;^^) bound comes from evaluating the continuous functions in the main 
and oscillatory terms att = i/rio. 

For Eg [Aipx/A(px'/] one uses the linear approximation of the angular shift to get 



A^e,x = Re[-(1 + zv)Ze,x] + 0{Z, 



I, a) 



and similarly for A'. After multiplying the two estimates and taking expectations, only 
the noise terms in a, ^^,a' contribute. Namely, with V = Vg, we have 

E, [Aifg^x^ipe^x'] = ^ E [(1 + 7]z)V + (1 + r]z)V] [{1 + r]z')V + (1 + r]z')V] + 0{k'") 
= ^ Re(l + zz')E\V\^ + ^ ReEV'^ + r] terms + 0{k'''^) 

where we used z = exp{i(fi^x) and z' = exp{iipi^x')- Formula (|69|) now follows from the 
asymptotics of E|V^p, EV"^. The last claim follows from the third moment asymptotics of 
X,Y. □ 



5.2 Continuum limit of the phase evolution 



The goal of this section is to show that the first stretch of the phase evolution converges in 
law to the solution of the SDE (f24)) . Typically, the phase evolves in an oscillatory manner, 
so we have to t ake advantage of averaging. Our main tool will be t he fo llowing propo- 
sition, based on 



Stroock and Varadhan ( 



197% and lEthier and Kurtz! (|l986h . which allows 



for averaging of the discrete evolutions. 



30 



Proposition 23. Fix T > 0, and for each n>l consider a Markov chain 

(X; G £= l...[nT\). 
Let Yl^{x) be distributed as the increment X^^^ — x given X" = x. We define 

b-{t, x) = n^lFL^ij (x)], a"(t, x) = nE[Y;:,,^ (x)^^:*] [xV]- 
Suppose that as n ^ oo we have 

|a"(t,x) -a"(t,i/)| + < c\x - y\ + o{l) 



n ( ^\ 1 31 



sup£;[|F/(x 



< cn 



-3/2 



(73) 
(74) 



xl 



and that there are functions a, h from M x [0,T] to respectively with bounded first and 

second derivatives so that 



sup 

x,t 



a^(s,x)ds— / a(s,x)ds 



+ sup 

x,t 



Vis.x) ds 



b{s, x) ds 



0. (75) 



Assume also that the initial conditions converge weakly: 



Then , < t < T) converges in law to the unique solution of the SDE 

dX = bdt + adB, X(0) = Xq. 

We will prove this in Appendix IA.3I The next lemma provides the averaging condi- 
tions for the above proposition. Recall that Atpi^x = ^\,e+i — ^£,\- 

Lemma 24. Fix A, A' and e > 0. Then for any £i < no(l — e) 



— ^ E [Aipi^xA^e,x' 



-- - V&A(t) + 0(/Xnno'/' + no"/') (76) 
' =y] = — a;, y) + 0(/i„no'''' + ng 



where t = i/riQ, the functions bx, a are defined in (fTOD , and the implicit constants in O depend 
only on e, f3, A. 

Proof. Summing (|68l) we get (|76l) with a preliminary error term 



— V Re(ei,, r/,) + — V Re(e2,, r/|) + 0{k;'^'] 

no Tin ^-^ 
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where the first two terms will be denoted Ci, C2- Here 

ei,i = {-v^{t) - ^g(^)/2)e-*^ e^,, = zg(t)e-2«/4 

and ki = uq — ii > cuq, where for this proof c denotes varying constants depending on e. 
Using the fact that v\, q and their first derivatives are continuous on [0, 1 — e] we get 

\ei/\ < c, \ei^£ - ei/+i| < cng ^ 

Thus by the oscillatory sum Lemma l37l 



ICil < c^{iJ,n/Vk + 1)^0^ + c(/i„/v^+ l)?^o ^ ^ c(/i„nQ^/'' + 



£=0 



Similarly, if we apply the same estimate for the second sum, we get 



^=0 



We could also estimate C2 by taking absolute value in each term. Using (|66l) together with 



(|47|l we get 
which leads to 

h-i 

1=0 

Using this bound for /i„ < 1 and the previous one for /in > 1 we get the desired estimate 
(|76l) . The asymptotics of the second sum follow similarly □ 

We are now ready to state and prove the continuum limit theorem. 

Theorem 25 (Continuum limit of the phase function). Suppose that riQ/n ^ \/{\ + v) with 
V G [0, 00]. Then the continuous function p(t) = p{nQ/n,tnQ/n) (see (1471 )) converges to a limit 
for which we use the same notation. Let B and Z be a real and a complex Brownian motion, and 
for each A G M consider the strong solution of 



A Rep' ^ Im(p2) 



2s Imp 2/3s2 



d^x = 

<^a(0) = TT. 

Then we have 

V^A.Lnotj ==^ Vxit), as n cx), 
where the convergence is in the sense of finite dimensional distributions for A and in path-space 
D[Q,l)fort. 
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Remark 26. From (|47|) we get that the Hmit of p{nQ/n, tuo/n) as no/n 1/(1 + z/) is 

V(l-t)(2z/ + l-t) 



Pit) 



V 



z/ + l- t z/+l-t 
with p{i) = \ \i V = oo. Thus equation (|77|) can be written as 

A 



Vl -t ci(^A = ^dt+J- Re{e-''^^dZt) 



V 



1 1 

^5 ^ 2/ Z/+ 1 -t 



'2(2z/ + 1 -t) 



t/E, (78) 



(5{v^\- t) 

where the last two terms are and 2(3'^^^dB, respectively when u = oo. 

Proof of Theorem^5\ It suffices to show that for any finite sequence (Ai, . . . , A^) and for any 
T < 1 the following holds on the time interval [0, T], 

We will use Proposition |23l For x G M'^ let 



be{x) 



X 



aAx 



UqE 



{A^){A^f\^=x 



Recall the estimates (|68|) and (|69|) . Since /i^/(4no) v, the functions 6a, « defined in (|70l) 
converge uniformly on [0, T] to h\,d which are also defined by (|70]) but in terms of the 
limit of p (recall that s is just \J\ — t). 
Using this with Lemma |24] we get that 



sup I / n^Yni)s\{p^)ds — I b{x,s)ds\ 0, 

xeR'',t<T Jo 

sup I / noainos\ix)ds — j d{x,s)ds\ 0, 

reR'i,t<T Jo 



(79) 



where 



a{x, t) ) = a(t, Xj, Xk) 



This means that condition (|75|) in Proposition |23] is satisfied. Because of ((TO]) and the mo- 
ment bounds we can see that (|73|) and (|74|) are also satisfied, thus {(px,inot\, ■ ■ ■ , ^\d,lnot\) 
converges weakly to the SDE corresponding to b{x,t),d{x,t). The only thing left is to 
identify the limiting SDE from the functions b{x, t), d{x, t). This follows easily, by observ- 
ing that if Z is a complex Gaussian with independent standard real and imaginary parts 
and iUi,uj2 E C then 



ERe(a;i-Z') Re{uj2Z) = E{uJiZ + ujiZ){uj2Z + 0O2Z)/4: = {uiUJ2 + oo2UJi)/2 = Re(ci;iCi;2). □ 
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Theorem |25l leads to the following corollary. 

Corollary 27. Let Wt he complex Brownian motion with standard real and imaginary parts and 
consider the strong solution of the following one-parameter family ofSDEs 

./T~i dax = \dt+ Re((e-^°^ - l)dWt). 

Then 

«L"oiJ,A =^ aA(t), as n ^ 00 (80) 

where the convergence is in the sense of finite dimensional distributions for X and in path-space 
D [0,1) /or t. 

Proof. If jj^/ (4no) converges to a finite or infinite value as n ^ 00 then the statement fol- 
lows immediately with Wt = e*'^"*^*^Z(. This implies that for any subsequence of n we 
can choose a further subsequence along which ([80)) holds, a characterization of conver- 
gence. □ 



5.3 Why are the limits in different windows the same? Universality 
and non-universality 

This subsection is meant to explain why the continuum limit of the relative phase function 
does not depend on the choice of the scaling window In order to do that, we will 
discuss a more general model where this is not necessarily true. 

The discussion of this section is not an integral part of the proof of the main theorem; 
the goal is to provide some additional insight for the results. 

A more general model. The following is a generalization of the model (|39l) . Consider 
random tridiagonal nxn matrices with diagonal elements Xq, Xi, X2, ■ ■ ■ and off-diagonal 
elements si, S2, S3, . . . and sq + Yq, si + Yi, S2 + Y2 . . . , see (|39l) . The random variables 
Xi, Yi are independent with mean approximately zero, variance approximately 2/P and 
a bounded fourth moment. The deterministic numbers si depend on n and are approxi- 
mately \pns{fln), where s{i) is a nonnegative, sufficiently smooth decreasing function on 
[0,1]. In the case of /3-ensembles we have s(t) = \J\ — t. 

We will try to understand the point process limit of the eigenvalues of these tridiago- 
nal matrices near /x„ where the scaling parameter /i„ will be in the interval [s(l)v^, s(0)v^. 
It turns out that in this more general setup the arguments of the previous two sections 
follow through essentially without change. This is the main reason why we expressed 
everything in terms of p, s, s, instead of using the sometimes more simple explicit values. 
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Figure 2: The definition of the scaling parameter ng 



We first have to identify the scaling around /i„ so that we have a nontrivial limit, for 
this we consider equation (|47|) . We define G [0, n] as the unique value for which 



(81) 



s{no/n) = Hn/V4:n. 



Then for i G [0, uq] the complex number G EI defined through (|47|) (see also (|48l) ) is of 
unit length and our scaling around /i„ will be given by (|50l) : 

A = /i„ + A/(2v/^). 

The subsequent computations, i.e. the introduction of the slowly varying phase function, 
the single-step asymptotics and the continuum limit of the phase evolution can be car- 
ried out the same way as we have done in Subsections 14.21 and 14.31 and Section |5l The 
assumption n^^jin^^ ^ as ri oo will ensure that the arising error terms are negligible. 

Thus, according to Theorem|25]if no/n converges to a constant as n — > oo then [noij , a — 
(fxi^t) where ^x{t) is the solution of ((77|). This gives the following limit for the relative 
phase function a^notiX- 



A 



sdax = -dt+ {{e 
Here s is given by the n — > oo limit of 



■tax 



l)dWt) 



«a(0) = 0. 



(82) 



Sn(t) 



' s{no/nty — s{no/ny 



no/n 



(83) 



There are various ways of interpreting equation (|82|) , perhaps the most intuitive is via the 
Brownian carousel which is already apparent in the discrete evolution, see Section |44l 
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The fact that the limiting equation depends on s through the function s explains 
two phenomena. First, in the /^-ensemble case s{t) = \J\ — t and thus (|83|) gives that 
s(t) = ^yl — t regardless of the value of the limit of riQ/n. This shows why all limits will 
be governed by the same stochastic differential equation, regardless on the choice of the 
scaling parameter. However, in the more general model, non-universality holds; the lim- 
iting stochastic differential equation ^2) depends not only on the limit of s but also on 
the scaling window. 

Second, consider a general s, with s'(0) < 0, and choose yu„ so that 

(2v^s(0) - iinW^ oo, and (2v^s(0) - iin)n'^'^ 0. (84) 

This means that scaling parameter yu„ is close, but not too close to the edge of the spectrum 
2s{Q)y/n. Since ^n/\f^ s(0),by (|8T]) we have no/n and 

no = s'{0)-^Vn{Vns{0) - /i„/2) + o{no). 

Thus tiq V"^^ ~^ 0/ so we can apply our previous results. From (|83l) we get 

s{t) = lims„(t) = cVT^, c = |(s2)'(0)|^/2 

which means that the limiting sde (|82l) is the same as in the /3-ensemble case, after a linear 
rescaling with a new parameter (3 = Pc^. This means that even for a general choice of s 



the point process limit of the eigenvalues in the scaling regime (|84|) is given universally 
by the Sine^? process. 

A similar statement of universality holds for a class of 1-dimensional discrete random 
Schrodinger operators with tridiagonal matrix representation. Consider the a symmetric 
tridiagonal matrix with diagonal and off-diagonal terms 

Xo, Xi, ... Zo/2 + so, Z,/2 + si, ... 

where Xi, Zi are independent random variables with mean approximately zero, vari- 
ance approximately and a bounded eighth moment. The deterministic numbers sg 
depend on n and are approximately y/ns{£/n), where s{t) is again a nonnegative, suf- 
ficiently smooth decreasing function on [0,1]. This gives the matrix representation of a 
1-dimensional discrete random Schrodinger operator. 

The analyze its spectrum, we first conjugate it with a diagonal matrix to transform it 
into a form similar to (|39|) . Choosing an appropriate diagonal matrix we can transform 
any off-diagonal pair {Zi/2 + s^,Z^/2 + sg) into (a^, {Yi/2 + sgf/ai) with any nonzero 
ai, while the diagonal elements stay the same, a simple computation shows that if we set 
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o-e = y/ns{i/n+l/{2n)) ~ -yns((^+ l)/n) then the off-diagonal entries above the diagonal 
will have mean approximately equal to y/ns{i/n), variance approximately equal to cr^ and 
a bounded fourth moment. Thus the previous results may be applied with P = 2/a'^. In 
particular near the edge of the spectrum (but not very near: see (|84)) ) the point process 
limit of the eigenvalues will be given universally by the Sine/j process. 

We would like to note that the point process limit near the edge of the spectrum 
(i.e. when — 2y^s{0))n^^^ converges to a finite constant) one gets the Airy^ process 
(see Ramirez, Rider, and Virag (|2007h . Section 5). This allows us to complete the proof in 
the general case with arguments analogous to the following section. To avoid excessive 
technicalities, we chose to focus on the beta ensemble case in this paper. We plan to treat 
the more general case in detail in a future work. 



6 Asymptotics in the uneventful stretch 

Section |5] describes the stochastic differential equation limit of the phase function on the 
first stretch [0, no(l—6)]. Here we show that in the limit, this stretch completely determines 
the eigenvalue behavior. 

6.1 The uneventful middle stretch 

The middle stretch is the discrete time interval [mi, 1712] with 

mi = [no(l -£)J, m2 = [n-fil/4-K{fifyi)\ (85) 

for e e (0, 1) and k > 0. The goal of this section is to prove that if a is close to an integer 
multiple of 2n after time mi then it changes little up to time m2- More precisely, we have 

Proposition 28. There exists a constant c = c(A, /?) so that with y = nQ^^{iin^ V 1) we have 

E[|a^,,A-a4,A| Aljj'^J < c(dist(a^,,A, 27rZ) + + y + ^-i), (86) 
for all K, > 0,e e (0, 1), A < |A| and mi < ii < £2 < m2. 

The first step is to estimate Aa^^A = 0!e+i,x — tt£,A using the angular shift Lemma [T6l with 
z = Zi^x defined in (|63|). For the finer asymptotics of the lemma, the condition \z\ < 1/3 
is needed. For this, we truncate the original random variables X^, Y^. For nii < £ < m2, 
introduce the random variables Xg, Yg which agree with X^, on the event 

\XA,\Y,\ < l7^s(£/no), (87) 
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and are zero otherwise; this event depends on n via s. By Markov's inequality and the 
fourth moment assumption (|40)) for Xg, Yi, this event has probability at least 1 — c(no — 



Summing this for £ < m2 shows that the total probability that the truncation has an effect 
is at most ck~^. This can be absorbed in the error term in (|86|) , so it suffices to prove 
Proposition l28l for the truncated random variables. 

To keep the notation under control, we will drop the tildes and instead modify the 
assumptions on X^, Y^. Namely, denoting k = uq — £ we assume the bounds (|87|) and the 
modified moment asymptotics 



moment 




2nd 


4th 






2/13 + 0{k'') 


0{1) 



which follow from the original ones (HO) and our choice of truncation. With p, q defined 
in (|67)) , this changes the moment asymptotics of Vi (|64]) the following way: 



(88) 



Proposition 29 (Single-step asymptotics for ae^x)- There exists k* = k*{i3, A) so that for every 
i < riQ — k* and \X\ < Xwe have the following. 



m 


ml 








4 


0{k-^) 


+ oik~^) 


^p(t) + o{k-') 


0{k~ 





— Re 

no 

+ — Re 

no 



'■«)r/ll +0{n-o'/'k-'/' + k-'/'ae,x) (89) 



0{k-'ae,x + no'^k'l') 
O{k''a,,x + no'k-') 
Oiae,xk-'). 



The functions vx = vx{i/no), q = q{i/no) are defined in 
between ag^x ^nd the set 2ttZ. 



(90) 
(91) 
(92) 

, and ae x denotes the distance 



Proof. By choosing a large enough k* > 1 we can assume that for i < n — k* 



X 1 
P - 10' 



which together with (|87|) guarantees that the random variable defined in (|63|) satisfies 
\Ze,x\ < Vs. The proof of the proposition relies on the evolution rule. Proposition I19|(ii)] 



Aai^x = ash((L 



-1, e"^'-'f]e) + ash(S,,o, e'^l^Vi, e'^'^'^Vi) + asli(S,,o, e^'^'-'f],, e'^^'^ve) 
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whose terms we denote C15 C2, Cs- First we show that Ci, C2 are small. By the definition ([551) 
of L we have 





1 A 




X 




uq 2s 







< 



10' 



This estimate with the third bound of Lemma [T6l gives 

Ci = ^l,-^e,x = 0{n-'^'k-'/'). (93) 
Applying again the third bound of Lemma [T6l with |^^,o| ^ and (|93l) gives 

For (^3 we use the first estimate of Lemma [T6l and note that in our case |i; — i^l equals 



|e*"AA _ i| < 



(94) 



Thus with Z = ZiQ we have 



(3 = - Re [{e-''^^'^ - e-''^'''')r]i{Z + zZV2)] + Re [i{e-^''^''^ - e-^^^^-")r7^ZV4] + C(d^,A^^). 

(95) 

Since Z is independent of JF^ and cee^x ^ ^£ / the error term becomes 0{ce£^xk ^^^) after tak- 
ing conditional expectation. The definition ((63]) of a and the moment bounds (|88l) imply 
that replacing EZ and EZ^ by vx{i/no) and q{i/no) gives error terms of order 0{k~'^). Be- 
cause of (|94|) and 

|g-2i<^,,^ _ g2ic^,,o| ^ |g2ia,,^ - 1| < 2d£,A (96) 

we get (|89l) . Using the explicit form of f a and g and (|94ll96l) again, we obtain (|90l) . The other 
estimates follow similarly from the first-order version of ((95]) and Proposition l22l □ 

The following lemma relies on the careful use of single-step bounds and oscillatory 
sum estimates. We postpone the proof till Section [6^ 

Lemma 30. Recall the definition of mi, 1712 from (|g5l ). There exist cq, ci defending on A, (3 so that 
with y = nQ^^{^n^ M I) we have 

3/2 



< < Cii^k + Unk + k '''Vnlfc>M2/4 j , 

whenever k > cq, |A| < A, and mi < £1 < £2 < m2. Here k = uq — £. 
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The last ingredient needed for the proof of Proposition |28] is the following determinis- 
tic Gronwall-type estimate. 

Lemma 31 (Gronwall estimate). Suppose that for positive numbers xe, b^, c, integers i\ < £2 
and £ = £1 + 1, £1 + 2, . . . , £2 we have 

i-\ 

xe<^ + c+J2bjXj. (97) 

Then 



Xi. 



< 2 (xfj + c) exp j 3 bj j . 



Proof. We can assume £1 = 0. Let ye = xe — xe^i/2, so that we have 

^e = yi + ^ + ^ + --. + ^ + ^. (98) 
Then and the positivity of Xq and bj gives 

,,<c + x,. + |:(,,(,, + l|i + l^t2 + ... + _|l_), (99) 



where s = 60 + • • • + ^^2-1- Taking positive parts in (|99)) , and then summation by parts 
yields 

yt <{c + xos) + Y,hyt (100) 

with bj = bj + bj^i/2 + ... + 6£2_i/2^2-i-i_ go that 1 < £3 < £2- We now multiply 

the inequality (flOOb by ^^(l + fo^+i) ■••(! + bti~i) and sum it for 1 < £ < ^3 — 1. We add 
(|100)) again with £ = £3. After cancellations, we get 

^3-1 ^ /^2-l \ 

yt^<{c + xqs) n (1 + &i) < (c + xqs) exp I ^ &i I < (c + xos)e^*. 
3=1 \ j=i J 



Applying this inequality for all the y terms in (|98|) with £ = £2 we get 

Xi^ <2{c + XQs)e^s + xq <2{xq + c)e"'. □ 

Proof of Proposition l2gl For this proof, let a = a^^^x, and define a^,a'^ e 2'kZ so that [a^, a^) 
is an interval of length 27r containing a. We condition on the cr-field JF^^, so in this proof E 
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denotes the corresponding conditional expectation. We also drop the index A from a. We 
will show that there exists cq so that if k > cq, then with the quantifiers of the proposition 

Eja^j-a^l < ci((a - a^) + + y), (101) 
E|a^2-a^| < ci((a^ - a) + Ve + y). (102) 

The claim of the proposition follows from this by an application of the triangle inequality 
to the stronger bound. The additional condition k > cq is treated via the error term 
Lemma l30l provides the bound 

1-2 

|Ea£ — a^l < (a — a^) + c(?/ + i/i) + Ea^_i/2 + 6jEdj. 

Note that a never goes below an integer multiple of I-k that it passes (Proposition |19|(iii)[ ), 
so (Xi > for all i > ii- This means that for i > £i we have < ae — and with 
Xi = E|q;£ — I we have the bound 

e-2 

xi < {a - a^) + c{y + y/e) + Xi_i/2 + '^bjXj. (103) 
According to Lemma |30] we can bound the sum of the coefficients be as 

£2-2 

^(^^2 + /^"'^2 + 1) < C' 

which means that (|101|) follows via the Gronwall-type estimate of Lemma |3T| 

Next, we consider the first time T > ii so that ax — > 0. Proposition |19|(ii) breaks 
one step of the evolution of a into two parts, from ai to a}_^-^ and from a}j^^ to ae+i- It 
shows that a can only pass an integer multiple of 2% in the first part. Since the first part 
is non-random, even the time T — 1 (and not just T) is a stopping time adapted to our 
filtration. The overshoot can be estimated in two steps. By ([93)) , and the fact that k > cq 
we have 

E [{a*T - a^)l(T < £2)] < cn~'^'. (104) 

By the expected increment bound (|90|) and the strong Markov property applied at T — 1 
we have 

E [{ar - «t)1(^ < ^2)] < crio'^'- (105) 

This gives 

E{ae,-a<>)+ = E [l{T < i2)E [a^, - a^l^r]] 

< Ci(E[(«r-a^)l(r<4)] +V^ + y) 

< c[iV^ + y), (106) 
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where the first inequality uses (|101[) and the strong Markov property, and the second uses 



(fT04l[T05ll . To prove (IT02ll first note that Lemma |30] also gives 



i-2 



|Ea^ - a9\<{a^ - a) + c{y + y/e) + Ed|._i/2 + ^ bjEaj. 

j=ii 

Then by (|106|) and the identity \a\ = —a + 2a^ we get 

E|a£-a^| < |Ea£-a^| +2E(a£-a^) + 

1-2 

< (a^ - a) + c(?/ + Ve) + Ed£_i/2 + ^ 6jEdj 



Since < \ae — a^|, the inequality (|103|) follows with = Eja^ — a^|, and the Gronwall- 



type estimate in Lemma |3T] implies (|102|) . □ 
6.2 Bounds for oscillations in the middle stretch 

This section presents the proof of Lemma |30l, isolated as the most technical ingredient of 
the proof in the previous section. We start with a bound on the mixed differences. 

Lemma 32. There exists an absolute constant c so that for i < n — k* (with k* as in Proposition 
|29D we have 

\Ei[Ae'^'-' - Ae'^'-^W < ck'^ae + cn~''^k~''^ 
and the same inequality holds with e^**^ replacing e**^ on the left-hand side. 

Proof. The left-hand side equals 

|e^^^,o [(gia^.A _ _ 1) + (e^^"^A _ i)(e^^'^Ao _ i) + (g^Aa,,, _ | 

< |e^"^-^ - l||E,[e^^'^^.^ - 1]| + E,|A«,,aA<^,,o| + |E,[e*^"^-- - 1]|. 

The statement now follows from ((68]), Proposition |29l the bounds ((94ll96l) and the bound 

|E[e*^ - 1]| < \E[e'^ -iX - 1]| + |E(iX)| < E\e'^ - iX - 1| + |EX| < EX^ + |EX|. 

The inequality involving e^*"^ can be proved the same way. □ 

Now we are ready to prove Lemma l30l 

Proof of Lemma l30l We will drop A in ai^x, and condition on the cr-field J^n^. Let E denote 
the conditional expectation with respect to this cx-field and let xi = Ed^ . We have 



l>2~l 



|EK-a,,]| < ^E[E(Aa,|j-,)] 



(107) 
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Let 



no no 4 



By the single-step asymptotics (|89l) the right-hand side of (|107[) can be bounded by 



^2"! r ^2-1 ^2-1 h-1 

\^Re{gi^eVe)\ + \^'Re{g2,e'nf)\ + \ ^ Re(5(2,(? 77^) | + c ^ /c '^''^x^ + c ^ no'''''A; 
i=ii i=ii 1=1" +1 i=ii i=ii 

with the usual notation k = no — i- We call the terms d, C2/ Cs/ C4/ Cs- Note that C2/ Cs 
come from a single sum cut in two parts at t = no — L/^n/4j , and one part may be empty 
Clearly, we have Cs < cy^, and ^4 is already in the desired form. By ((94|) and the bounds 
(|65l[67l) on g, v we have 

c c 

\giA < — I^^A + ^g/2|x£ < 

77-0 K 



Lemma |32] with t+ = (£ + l)/no gives 



no 

< ck~'^xe + cno^^^k 



where we used the notation A^/ = /((£ + l)/no) — f{i/no). The oscillatory sum Lemma 
|37| gives 

£2-2 

Ci < c{fin{no - i-iY'l' + {no - £2)"')x£,_i + c 5^(a;,A;-2 + no'/'A;"/')(/inA^"/' + 1) 

£=£1 

where the coefficient Ye is achieved by choosing a large enough cq. We continue 

£• r £* 

C , , . , , -3/2 



C2 < X^|5'2/| < — klX^a;;? < c xi fink 



The term ^3 is handled by Lemma |37l with gj = g2j. Standard bounds on q, q' and Lemma 
give 

\ge\ < ck'^xi, \gi - gi+i\ < ck''^Xi + cnQ^^^k 
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hence from Lemma |37| we get 



Cs < 



i=i*+i 

£2-2 



< 



i=t*+i 



if Co is chosen sufficiently large. The claim follows. 



□ 



6.3 Why does the right boundary condition disappear? 

The goal of this section is to show that the phase evolution picks up sufficient randomness 
that will neutralize the right boundary condition of the discrete equations. 

Proposition 33. Letm = [n— /i^/4— K(/in^Vl)J and suppose that k 00 with nQ^K{j/n^\/l) — 
0. Then ip.m,o modulo 27r converges in distribution to Uniform{0, 27r). 

Proof. We will show that given e > 0, every subsequence of indices has a further subse- 
quence along which Lpmfi modulo 27r is eventually e-close to uniform distribution. So we 
first pick an integer r = t{e) and show that along a suitable subsequence, the conditional 
distribution given Tm-Tt, of v^m.o — V^m-T^,o converges to Gaussian with variance tending 
to 00 with r. Here the scaling factor is ^ = \fi{^n V 1)J . Since a constant plus a Gaussian 
with large variance is close to uniform modulo 27r, the claim follows if we let r go to 00. 

To show the distributional convergence, we apply the SDE limit Theorem |25] to the 
evolution of from time m — on. To adapt to the setup of the theorem we introduce 
the new parameters 

fi = n-m + T^, /ifi = /i„, rio = n - /i?/4 - 1/2, (pe,\ = (p£+m^ri,\- 

By assumption, we have nQji^^ —>■ 00. We pass to a subsequence so that no/h has a 
limit 1/(1 + z/) G [0, cxd], so Theorem |25l (trivially modified to allow general initial con- 
ditions) applies. The result is that (fitno\,o has an SDE limit given by (|78l) with A = 0. 
Thus (pm,o — ^m~T^,o converges to a normal random variable which does not depend on 
the initial value <^(0). Its variance is given by integrating the sum of the squares of the 
independent diffusion coefficients on the corresponding scaled time interval: 




l-(l+r) 



-1 




dt > -log(r + 1), 



which goes to 00 with r, as required. 



□ 
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6.4 The uneventful ending 

This section is about the last part of the recursion, from 

m2=[n- 111/ A - K{iif V 1)J 

to n where k > is a constant. The goal is to show that nothing interesting happens on 
this stretch. More precisely, we show 

Lemma 34. For every A G M and k > as n —>■ oo we have t^?® ^ ^ — ip'^^ q ^ m probability. 

Fix K and A. We will show the convergence by showing that any subsequence has a 
further subsequence with the desired limit. Because of this, we may assume that the limit 
of fin exists. We will consider two cases: lim yU„ < oo and lim yU„ = oo. 

Proof of Lemma\3^in the case when lim /i„ is finite. 

In this case we can assume that n — m2 is eventually equal to some integer ^. Also, 
converges to a unit complex number p with Im p > 0. By (|45|) we have 

^ti,x*Qn-^-i = o*R;\ • ■ ■ R;!^, (108) 

where R^l^ x = ^n-j^n-j aQ(^)^- Consider the components of the product on the right- 
hand side of (|108|) . The elements Ln-j,A are deterministic (see ^T} ) and as functions on 



i' - the lifted unit circle - they converge uniformly to non-degenerate limits that do 



not depend on A. (Here we also used s„_j = y j — V2.) In the same sense, we also 



have T„_^ A(Im(p)^^, — Rep). Because of the moment bounds (|40)) we may find a 
subsequence along which X„_i, . . . , and Yn-i, . . . , Yn^^ all converge. Then (using 
the definition (|43|)) it follows that the random elements W„_j converge as functions for 

Since all of these limits are non-degenerate and the dependence on A disappears, we 
have 

Remark 35. For the second case, we review some of the results of Ramirez, Rider, and 



Virag (|2007l) , henceforth denoted RRV, about the eigenvalues of the stochastic Airy oper- 
ator. The paper considers the eigenvalue process A„ of the random matrix M (see (|2T|)) 
under the edge scaling n^/''(A„ — 2y/n). By Theorem 1.1 of RRV, the limit is a point pro- 
cess S given by the eigenvalues of the stochastic Airy operator, the random Schrodinger 
operator 
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on the positive half-line. Here h' is white noise and the initial condition for the eigenfunc- 
tion / is /(O) = 0, /'(O) = 1. By RRV, Proposition 3.5 and the discussion preceding RRV, 
Proposition 3.7, 

S is a.s. simple, and for every x G R, we have P(x G S) = 0. (109) 

The proof is based on the observation that after appropriate rescaling the matrix M acts on 
vectors as a discrete approximation of Hp. The initial condition /(O) = 0, /'(O) = 1 comes 
from the fact that the discrete eigenvalue equation for an eigenvalue A = + n~^/^v is 
equivalent to a three-term recursion for the vector entries wi^y (c.f. (|4T|) and Remark |20]) 
with the initial condition w^^y = and Wi^y ^ 0. 

By RRV, Remark 3.8, the results of RRV extend to solutions of the same three-term 
recursion with more general initial conditions. We say that a value of v is an eigenvalue 
for a family of recursions parameterized by u if the corresponding recursion reaches in 
its last step. Suppose that for given C, G [— oo, oo] the initial condition for the three-term 
recursion equation satisfies 

where p„ := Wi^y/wQ^i, does not depend on v. Here the factor n/''' is the spatial scaling 
for the problem (RRV, Section 5). Then the eigenvalues of this family of recursions con- 
verge to those of the stochastic Airy operator with initial condition /(0)//'(0) = C,. The 
corresponding point process will also satisfy (|109|) , see RRV, Remark 3.8. 



Now we are ready to complete the proof of Lemma |34l 
Proof of Lemma |3H in the case when lim /x„ = oo. 

Without loss of generality we assume A > 0. Fix a G M' and let B denote the event that 

3^*Qm2-i ^ ^ mod 27r for x G V^m2,o 
It suffices to show that P(i?) 1. Indeed, by considering a subdivision of the unit circle 



into arcs of length at most e at points e \ if the event (|110[) holds for each 6j then 



|v^m2,A*Qm2-l V^m2,0*Qm2-l| ~ |v^m2,A V^m2,o| 

cannot be greater than e. Taking e ^ completes the proof. 

Equation (|54|) translates B to an event about (pf^^. More specifically, by Proposition [TTI 



B is the event that the one-parameter family of recursions parameterized by v 

(pt+l,v = ^e,u*'Ri,iy, i>m2 
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with initial condition 

0m„u = e.T-^\ (111) 

does not have an eigenvalue in the interval [0, A]. This recursion is determined by the 
bottom right n2 x n2 submatrix of M(n)^^^^ (|39l) , where n2 = n — m2. Thus the recursion 
is in fact the discrete eigenvalue equation for M(n2)^*^"^^ with a generalized initial con- 
dition. This can be transformed back to the discrete eigenvalue equation for M{n2) with 
the corresponding initial condition. Let u = U^(e*^) G R be the point corresponding to 
9 E M'. Then (|111|) translates to the initial condition 



^m2 J ) 



for the eigenvalue equation of (see dHJ) and by Remark|20]to the initial condition 

X{n-m2-l)l3 X(n-m2-l)f3 /t f \ \ -o f W /i n\ 

Pm2,v = rm2M —7^ = (Im(p^jM + Re(p^J) (112) 

VP 5^2+1 \/P\Jl -1712- VV 

for the eigenvalue equation of M. As /^^ oo, we have 

n2 ^ oo, fil^l A = n2 — Kn^^ + o{nl^), and = 1 + ^V^^2^''^ + '^(^2^''^)- (113) 
Since xi converges to 1 in probability as £ — oo, (|112|) and (|113|) imply 



'^2^''' {Pm2,u - 1) ^ ^^^U ^ =: C- 

This means that the limit of P(-B) can be related to the limit point process S^. The interval 
[0, A] corresponds to 2p„ + [0, in our scaling ((50|). In the edge scaling correspond- 

ing to n2, the length of the remaining stretch, this turns into the interval 

(2p„ - 2nf ) + [0, Ano'/72] ^ -k + [0, 0] 

where the convergence follows from (|113|) . 

For 5 > let Bs be the event that the discrete eigenvalue equation for M{n2) with 
initial condition (|112|) does not have an eigenvalue in the interval 



2nl^ + n2^^^ {—K — 6, —k + 6). 
By Remark |35l, for any fixed S we have 

lim sup P(-B) < lim sup F{Bs) < P(5^ doesn't have a point in — 5, k + 5]). 



Since this holds for all 6 > 0, the fact (|109|) gives limP(i?) = 1, as required. □ 
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A Tools 



A.l Angular shift bounds 

The objective of this section is to prove Lemma [161, which relies on FactlTSl 

Proof of Fact^5\ The general form of such a transformation is given by WcT = e^°'{w — 
cr)/(l — (Tw), where a is the pre-image of 0. We may assume a = since post-composing 
T with a rotation does not change the quantities in question. Using the definition of 
ash(T, f , w) and the fact that |w| = = 1 we have 

ash(T,t;,t/;) = Argj^^.) [j^-—- ) " ^^S[o,2.)i^M- 



The additivity of Arg mod 27r proves (|34|) mod 2n. By definition, ash is continuous in T 
and so also in a. Since \a\ < 1, the right-hand side of (|34|) is continuous in a. As equality 
holds for (7 = 0, the proof is complete. □ 

Proof of Lemma^^ Recall that r.U = (i — r)/(i + r) maps the upper half plane to the unit 
disk, sending z to 0. By Fact [15] we have 

, N . f \ — iii + z) .\J) w\ . , , ziw — v) 

ash(T,t;,u;) = 2Arg ^— ^ = 2Arg (1 + x) , x ^ ' 



l~ + z).\])v J " 2i + z{l + v) 

If kl < Vs then we have \x\ < 1/2 so we can write asli(T, t>, w) = Re hy^^{z) with 



KAA = - log 1 + . =[w-v)\-z z + r/„,^ 

I \ 2i + z{l + v) J \ 4 

Here we use the standard branch of the logarithm defined outside the negative real axis. 

The second equality is Taylor expansion in z. To bound the error term, we write 

, /,/ / X ^ [w -v)p{z,v,w) 

nv,u,KZ} + z{l + v)f{2i + z{l + w)f ' 

where p is some (explicitly computable) polynomial, so the Taylor error term satisfies 

\riv,w[z)\ < —7- sup < c\z\'. 

<j! |2|<1/3>| = |«>|=1 \W — V\ 

This proves the quadratic approximation of the angular shift for \z\ < 1/3, and the other 
two estimates of (|37l) follow easily. 

For the case \z\ > 1/3, v = —1 we use the fact that |Arg(l + x)| < 7r|a;| to conclude that 

z{w — v) 



ash(T,f,iu)| < An 



i + z{l + v)/2 



4n \z{w -v)\< 47r3'^"^ - v) 



for any d > 1. Using \z\ > 1/3 we get that the main terms on the right-hand side of (|37|) may 



also be bounded by q|w — f | l-^l'^ and from this we get upper bounds of (|38|) as well. □ 
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A.2 Oscillatory sums 



Recall from the definition (|53|) that T]i is a unit complex number with a rapidly oscillating 
angle. Lemma |37]below will show that this oscillation has an averaging effect in sums. In 
order to prove that we first need the following harmonic analysis lemma. 

Lemma 36. Suppose that 2n > 9i > 62 > . . . > Om > and let Si = Yl]=i ^j- Then 



max \J2e-^\<c{9-^ + i27r-9,)-'). 



Proof. We first consider the case when 61 < n. Using second order interpolation we can 
construct a differentiable function s{x) on [l,m] with s(£) = si for 1 < i < m for which 
the derivative s'{x) is monotone decreasing derivative with — tt < s'(x ) < —9^ /2. 



Our proof is base d on the following lemmas of van der Corput (see Hilld (|l929l) for the 



first and 



SteinI (|l993h . Chapter VIII, Proposition 2 for the second): 



(i) If s{x) has a monotone derivative with < tt for x E [a, b] (with a,b E Z) then 
the difference of X]^=a ^^'^ la ^^^'^^^dx is at most 3. 

(ii) If s'{x) is monotone and > p on an interval [a, h] then | f^e'^^^^^dx\ < 3p^^. 

Since for our function tt > \s'{x)\ > 9.^/2 for x E [1, m] we may apply these lemmas to get 
the bound | Y!j=i e*'' I < c^„^ 

Consider now the case 2tx > 9i > ti . Let be the largest index with 9i* > tt, then 



£ 



(114) 



j=i u=i j=i u=i j=e*+i u=e*+i 

The second sum can be bounded by c 6^"^ using the first half of our proof. To bound the 
first sum, note that 

^ exp ^ 6'„] = ^ exp ' Ou = 2n - 9^+1^^ 

j=l u=l j=l u=l 

and for i < £* we have 

71 > 271 - 9i, > 9i > 92 >...> 9i > 271 - 9i > 0. 
Thus the first half of the proof can be applied again to get the bound c(27r — 9i)~^. □ 

The following lemma describes the averaging effects of the oscillating unit complex 
numbers 77^. 



49 



Lemma 37. Let E Cfor i eN and io < ii < uq. Then 



h ei-i 



1 - ge\ 



(We used the shorthanded notation k = no — £, ki = no — ii and ko = no — io-) 

Proof. For d = 1,2 we introduce F^j = X]m=£o with -Fd/o-i = 0. By partial summation 



Yl 3jVj = Fd/,gi, + Y ^ddigj - gj+i)- 

From (|49]) we get the estimates 

Axgpi < /i^^A;'^', and ti/2 - Argp^ < /i„A;"''''. 
Together with (|53|) this means that we can use Lemma |36] to get 

\Fi^t\ < c (p-nki^'^ + \F2^t\ < c (p-nki^'^ + /i^^^A:, 

This with (|115|) implies the lemma. 



(115) 



□ 



A.3 A limit theorem for random difference equations 

Proof of Proposition |23l Let || ■ ||oo denote supremum norm on [0,T]. For a two-parameter 
function / and x E M let J denote the integral J = fl f{s,x) ds. We recycle this 
notation for a function X : [0, T] ^ M to write T f^xif) = Jq f( s,X(s))ds. 



The proof of this proposition is based on Theorem 7.4.1 of lEthier and Kurta (|l986l) , as 
well as Corollary 7.4.2 and its proof. (See also lstroock and Varadhan ( 1979 ).) It states that 
if the limiting SDE has unique distribution (i.e. the martingale problem is well-posed) 
and also 



1b",X" ~ ^ b.X^^Woo ^ 0, 

II P 



IX 



-I 



a,X" oo 



0, 



for every e>0 supP(|r/(a;)| > e) — ^ 0, 



(116) 



(117) 



then X" =^ X. The theorem there only deals with the case of time-independent coeffi- 
cients, but adding time as an extra coordinate extends the results to the general case. 
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Because of our assumptio ns on a and b the well-po sedness of the martingale problem 



follows from Theorem 5.3.7 of Ethier and Kurta (|l986l) (see especially the remarks follow- 
ing the proof), and even pathwise uniqueness holds. Condititon (|117|| follows from the 



uniform third absolute moment bounds (|74)) and Markov's inequality. Thus we only need 
to show (|116|) as well as the analogous statement for a, for which the proof is identical. 
We do this by bounding the successive uniform-norm distances between 



I, 



X 



X 



where X'l'^ = X^^^/^j with K = \nT/L\, and = X"^fj. In words, we divide 

[0, \ nTW into L roughly equal intervals and then set X"'^ to be constant on each interval 
and equal to the first value of X"^ occurring there. 

If a function / takes countably many values fi, then for any h we have 



IX 



h,f\\oo 



< 



Since X"'^ takes at most L values, we have 

\\Zb",X"'^ — Xft^X^.-f- llcxD = llcxD ^ -f' sup ||X5n_^^^.||oo = Lo(l) 

X 

by (|75|) where o(l) is uniform in L and refers to n — > oo. From (|73|) , the other terms satisfy 



X J5^7i.L — 1 l^n 



Hoc 



< T||6"(-,X"'^(-)) -&"(-, X"( 



< cT||X"-X 



n,L I 



The same holds with b replacing 6". It now suffices to show that 



EllX"'^ -X"| 



Esup|X;'^-X,"| < f{L) 



(118) 



uniformly in n where f{L) — as L — > oo. The left-hand side of (|118|) is bounded by 



£—1 £ 

Esup|x7-i &n(x;) -x;-^! + Esup|^ ^ 6(x 



and the second quantity is bounded by T sup^^ 1 6" (x) | /L. The first quantity can be written 
as EM* where 



M* 



max M* 

i=0,...,L-l 



M* 



max I Mj 



=0,...,K-1 



Mi/ — XiK+t — XiK 



1 



iK+k) 



k=0 
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Note that for each i, Mi^£ is a martingale. For any martingale with Mq = we have 



Emax|M,f < ceI VE[(Mfc - Mfc_i)2|^,,_i] 

k<n I ^ 



k<n 



3/2 



< cn^/^maxE[|Mfc - Mk-i\^\J='k-i] 

k<n 



The first step is the Burkholder-Davis-Gundy inequality (see iKallenbergI (|2002h . Theorem 
26.12) and the second step follows from Jensen's inequality. Therefore (|74)) implies 

E[|M;| V^l] < c{n/Lf'n''/' = cL-''\ 

which gives the desired conclusion 



L-l 



(EM 



< 



j=0 



Letting first n ^ oo and then L ^ oo gives (|118|) and (|116|| . 



□ 
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